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We present a method for transferring bosonic atoms residing on the lowest s-band of an optical 
lattice to the first excited p-bands. Our idea hinges on resonant tunneling between adjacent sites 
of accelerated lattices. The acceleration effectively shifts the quasi-bound energies on each site such 
that the system can be cast into a Wannier-Stark ladder problem. By adjusting the acceleration 
constant, a situation of resonant tunneling between the s- and p-bands is achievable. Within a 
mean-field model, considering 87 Rb atoms, we demonstrate population transfer from the s- to the 
p-bands with around 95 % efficiency. Nonlinear effects deriving from atom-atom interactions, as 
well as coupling of the quasi bound Wannier-Stark states to the continuum, are considered. 

PACS numbers: 03.75.Lm,03.75.Nt,05.30.Jp,67.85.Hj 



I. INTRODUCTION 

With the recent experimental progress, systems of ul- 
tracold atoms in optical lattices have become a princi- 
ple playground for the study of correlated many-body 
systems typically appearing in condensed matter theo- 
ries [l|. For example, the ability to cool down atomic 
gases to very low temperatures, having access to the po- 
tential parameters as well as to the atom-atom interac- 
tion strength, made it possible to detect, for the first time 
experimentally, the Mott-superfluid phase transition in 
these systems Q- 

To date, most of the experiments have been confined 
to the lowest energy band of the optical lattice, the so- 
called s-band, and effects arising due to higher bands 
could be safely ignored. Contrary to the s-band, the on- 
site wave-functions, or Wannier functions, on the first 
excited bands, the p-bands, have a single node in some 
spatial direction giving them an orbital character. In 
the direction of the node, the onsite functions are more 
extended and therefore the tunneling coefficients in this 
direction is greater than in the perpendicular directions. 
Due to this anisotropic structure, the physics on excited 
bands allow for new phenomena not present on the lowest 
band, like supersolids or other types of novel phases Q- 
Moreover, the p-bands are two-fold or three-fold degen- 
erate (in two or three dimensions respectively), and con- 
sequently enable non-trivial phases even within the insu- 
lating phases Q. 

In fermionic systems, the p-bands can be excited sim- 
ply by considering filling factors larger than unity such 
that the lower s-band is completely filled 5]. For ultra- 
cold bosons, which we consider in this work, populating 
the p-bands implies an excited non-equilibrium state and 
the system is thereby energetically not stable. Even so, it 
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turns out that the life-time for bosons on the p-bands is 
long in comparison to characteristic time-scales for tun- 
neling between neighboring lattice sites 0, Q. Thus, de- 
spite decay, non-trivial physics can indeed be experimen- 
tally realized on the p-bands, for example condensation 
into non-zero momentum states has already been seen 0- 
9] . The experimental methods used for preparing bosons 
on the p-bands include; Raman transitions via a third 
band [10] , adiabatically ramping up a constantly moving 
lattice iQ, or sudden switching of superlattices @Jp]. In 
the last scheme, loading of both p-band bosons [8[ and 
/-band bosons [9| were achieved. 

In the present work we analyze, on a mean-field 
level, an idea somewhat similar to the one employed in 
Refs. [a, |9j, but for a regular monochromatic lattice, i.e. 
a lattice with only one characteristic wave-length A. We 
suggest a sudden tilting of the lattice via acceleration of 
it [3, [HI [Hi • The dynamics in the tilted system is con- 
veniently studied in the Wannier-Stark picture where we 
form a basis of localized states at each lattice site. The 
corresponding energies of these states are the Wannier- 
Stark energies, which form resonances in the full con- 
tinuous spectrum of the system. In the regime of deep 
lattices, we impose a tight-binding approximation where 
transitions appear only between neighboring lattice sites. 
By adjusting the tilting, a resonance coupling between 
the s- and the p-bands occurs. Suddenly switching off the 
acceleration/tilting after half a Rabi period leaves most 
of the initial s-band atoms in the p-bands. Transitions 
into higher excited bands are suppressed due to the off- 
resonant transitions with these bands. Including the s-, 
p-, and d-bands, we consider the time-dependent Hamil- 
tonian which takes the switch on/off into account. We 
analyze both the one and the two dimensional cases. It 
is furthermore shown that interaction between the atoms 
leads to a self-trapping effect and hence a decrease in 
the loading efficiency. This, however, becomes impor- 
tant only for very strongly interacting systems and should 
therefore not be a problem as long as one considers a 
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weakly interacting gas. Throughout we use experimen- 
tally relevant parameters and demonstrate a transfer of 
atoms from the s- to the p-bands of more than 95 %. 



II. MODEL SYSTEM 
A. Multi-band Wannier-Stark ladders 

The spectrum of a Schrodingcr operator on the form 

(1) 



h 2 d 2 



H = -~ + V L ^(kx) 



has a band-structure E v (q) characterized by a discrete 
band index v = 1, 2, 3, ... and a quasi momentum 
—2%/X < q < 2ir/\. Here, Vl is the potential depth, 
m will be the mass or our atoms, and k = 2ix j\ the wave 
number where A is the lattice wave-length. For the lat- 
tices that we consider in the present paper, we label the 
first three bands s, p, and d respectively [IJ- For higher 
dimensions there are multiple p- and <i-bands. When- 
ever the depth of the lattice potential is increased, the 
band- widths \E V {Q) — E v {2tt / X)\ decrease. In the limit of 
infinite depth, the bands become completely flat and rela- 
tive to the onsite oscillator energy approximately equidis- 
tant, mimicking the spectrum of a harmonic oscillator 
E v = (47r/A) J Br dq E u (q) oc v. In such a situation, 
particles residing in the lattice are immobile. 

For a constant force with strength Fq applied to our 
system, 



2m dx 2 



+ Vl sin 2 (fcir) — F x, 



(2) 



it is practical to define the Wannier-Stark energies [l4| 
which approximate 



E v +jF ^-iT 



(3) 



with j being an integer labeling the sites and gives 
the decay of the z/th band into the continuum of states. 
For each z/, the energies form a so called Wannier-Stark 
ladder. Given j, the low lying energies correspond to 
quasi bound states localized at site j. A schematic pic- 
ture of the Wannier-Stark energies is presented in Fig. [1] 
showing the first three bands/ladders. In this example, 
we have taken Vl = 20E r , where E r — h 2 k 2 /2m is the 
recoil energy. The linear force is taken such that for the 
real parts E2j = i-e 



2(E 2 - gi) 
= - x c, 



(4) 



where c is a dimensionless parameter that we have in- 
troduced for later purposes. More precisely, in the full 
system an infinite number of states are coupled while only 
a fraction of them couple resonantly. Population transfer 
occurs predominantly among these resonant states and 



off-resonant states can be, in a first approximation, ig- 
nored. Nevertheless, the virtual processes related to cou- 
pling to far detuned states induce small energy shifts in 
the states of interest, and consequently shift their mutual 
resonance conditions. An expression for these shifts can, 
in principle, be obtained by adiabatic elimination of off- 
resonant states, but this is not manageable for the system 
sizes we study, and instead we introduce the constant c 
in (j4]). In other words, when necessary c adjusts for the 
virtual energy shifts originating from off-resonant pro- 
cesses. In the figure c = 1. For these parameters, the s- 
and p-bands are almost completely flat while the d-band 
has a small width. Even though the bands (especially the 
e?-band) have a small width, neighboring p- and cZ-bands 
are mutually off resonant for any quasi momentum q. 




FIG. 1: A tilted periodic sine potential with its correspond- 
ing Wannier-Stark energies (here we have not averaged the 
energies over the Brillouin zone in order to demonstrate the 
band widths). The potential depth is Vl = 20E r . 



The imaginary part T v of the Wannier-Stark energies 
can be made relatively small for the s- and p-bands as- 
suming sufficiently deep lattices, as for example in Fig.[T] 
When the coherent coupling between the two resonant 
states are large compared to the decay rate, the two 
states can perform numerous of Rabi oscillations before 
they decay. This is the present idea for loading of bosons 
into the first excited bands. The energies E v can be 
obtained by direct diagonalization of the untilted Hami- 
tonian. The imaginary parts T v can either be found by 
diagonalization of the full evolution operator [14| or es- 
timated using Landau-Zenner theory [3, EH- For our 
set of parameters, we found it numerically difficult to 
get r„ utilizing these methods. Instead we estimate T v 
for the two lowest s- and p-bands by performing a wave- 
packet propagation of the full Hamiltonian and calculate 
the decay of population out from the first two Brillouin 
zones 



B. Tight-binding model for accelerated lattices 

As a basis for our analysis, we will consider 87 Rb atoms 
with mass m (= 1.44 x 10 -25 kg) trapped in an isotropic 
harmonic trap with frequency uj = 40 x 2-7T Hz. The atoms 
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are exposed to an optical lattice, with the wavelengths 
A = 795 nm and the depths Vl — 20E r the same in all 
spatial directions [13] • These are typical experimental pa- 
rameters, even though the lattice depth is relatively deep, 
however still well within experimental reach. We have 
found numerically that a deep lattice is preferable for 
our loading scheme. Only the one- and two dimensional 
cases will be considered in this work (the generalization 
to three dimensions is straightforward, but numerically 
demanding). The derivation of the Hamiltonian will be 
for the two dimensional case. The lattice is formed by 
counter propagating waves in such a way that the rel- 
ative frequency Acjj,(i) of the two waves can be tuned 
in time t. Such lattices, with Aw(t) = k [at 2 — bt) (a 
and b are constants), are typically used for the study of 
Bloch oscillations of cold atoms [ll|. Within the mean- 
field regime, and assuming s-wave scattering between the 
atoms, the equation of motion for the order parameter, 
after a proper unitary transformation [l2j], reads 



_ d . , . 
ih—ip[x,t) = 



V 2 + V tr (x) + V lat (x) + Fx 

Am 



Un 



Wx,i)| 5 



^(x,i), 



( 5 ) 

where F = ma is the amount of tilting of the lattice, and 
Viat(x) = V L [cos 2 (fcc) + cos 2 (fcy)]. 

Initially we will assume F = and the atoms to be 
prepared in the lowest s-band. The acceleration of the 
lattice is switched on during a time T s and stays on for a 
time T r and then switched off during T s again. To mimic- 
such an experimental ramping, for our numerics we take 
the force as 



tanh 



t - 0.1T r 



tanh 



t - l.lT r 



(6) 

with Fq given in Eq. (j4j. This F(t) has a smoothened 
square pulse shape, starting at t = 0.1T r and ending at 
t = l.lT r , and with a turn on/off time T s . As explained 
above, the parameter Fq and the time T r are chosen such 
that the lowest Wannier-Stark energy is resonant with 
the first excited states of a neighboring site, and the sys- 
tem undergoes half a Rabi cycle. Due to the anharmonic- 
ity of the lattice potential, the second excited Wannicr- 
Stark energies will be off-resonant with neighboring first 
excited Wannier-Stark energies. Nevertheless, for small 
detunings, higher bands may still be populated and in 
our calculations wc therefore include all the s-, p-, and ri- 
bands. In two dimensions, the first excited bands are two- 
fold degenerate and the second excited bands three-fold 
degenerate. We will specifically pick the lattice depth to 
be 20 recoil energies, Vl = 20E r . In this case the low ly- 
ing bands are relatively flat, but still the anaharmonicity 
is present. This was demonstrated in Fig. [TJ 

Apart from restricting the numerics to contain only 
the s-, p- 1 and ri-bands, we will further impose the tight- 
binding approximation. This means that we consider 



hopping only between neighboring sites. Both these ap- 
proximations should be justified for relatively large lat- 
tice depths Vl- As the localized Wannier-Stark states 
at each lattice site i = (i x , iy) we take the numerically 
obtained Wannier states This is a good approxi- 

mation for the low lying states of deep lattices [1J|. In 
order to derive mean-field equations of motion in the re- 
stricted Hilbert space, we expand the condensate order 
parameter as 



(7) 



where a = 0, x, y, xy, xO, Oy and i runs over all lat- 
tice sites. The a-index represents six possible Wan- 
nier states (s-band: "0", p-band: "x" and "y", and ri- 
band: "xy" , "Ox" , and "Oy" ) . Normalization is such that 
No = j rix |i/>(x, t)\ 2 , where A*o is the total number of 
atoms. The energy expressed in the coefficients ipf(t) 
becomes, 



a\2 



>/3 (SJ) 



(8) 



a, /3, 7,<5 i 



where 



rix wf (x) 



2m 



V 2 +y /Qt (x)+Fx 



wf (x), 



U a/ 3 7 s = U scat J dxwi(x)w^(x)w^(x)wi(x). 

(9) 

The summation runs over nearest neighbors, we 

have assumed the atom-atom interaction to be non-zero 
only within one and the same site, and asterix * denotes 
complex conjugation. The adjustable parameter U sca t 
is determined by the system dimensions and the s-wave 
scattering length a s . The effect of the trapping poten- 
tial is seen only in the shift of the onsite energies. This 
should be justified since the slope of V tr (x) is typically 
a thousand times smaller than \Fq\ within the atomic 
cloud and should therefor not induce any hopping be- 
tween sites. We further note that the lattice Hamiltonian 
Hint = p 2 /2m + Viat(x) only couples Wannier functions 
within the same band, i.e. {wf\Hi at \wf) — whenever 
a ^= (3, and as a consequence it is the tilting alone that 
induces the transition from the s- to the p-bands. 



III. LOADING IN A ONE DIMENSIONAL 
SYSTEM 

The expression for the energy (jSJ) is general for any 
dimensions. In one dimension, the excited bands are not 
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degenerate and we can label them by a = s, p, d. As we 
will demonstrate, the underlying physics is easily under- 
stood in one dimension. Nevertheless, some of the effects 
are more significant in two compared to one dimensions 
and vice versa. 



A. Two-site problem and the effect of interactions 
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FIG. 2: Populations P StP (t) in the s- (pink curves) and p-band 
(black curves) . In the upper two plots we compare the results 
obtained from the two-state effective model, Eq. ()12|) (dotted 
curves), with those obtained from the two-site six-state model 
defined in Eq. (|f 0|l (solid curves). In (a) U s = U p — and 
in (b) UsNq ss 4.3t$ and U P N « 3.1<*|. The effect of the 
interaction is manifested as self-trapping. The lower plot (c) 
displays the full iV-site model (N = 41) including a trapping 
potential. The interaction is set to zero in this case. In all 
examples Vl = 20E r , the tilt parameter c = I, and for the 
trap frequency in (c) w = 40 x 2n Hz. 



Effects deriving from atom-atom interactions, as well 
as the dynamics in general, can be extracted by consid- 
ering the most simple situation consisting of a lattices 
with only two sites similar to a double- well system. Here 
we consider the static model where the force F(i) = Fa 
is time-independent such that we expect Rabi evolution. 
Thus, we have a set of six coupled equations for the coefH- 
cients ^f(t), <fi(t), yi{t), f 2 (t), yl(t), V d 2 (t). The atom- 
atom interaction-free Hamiltonian H 2 reads 
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(10) 



where we have introduced the detunings Ay = £i,i — S2,j- 
Only the tp\ (t) and the <p^(t) states are degenerate. More- 
over, is typically 2-4 orders of magnitude smaller 
than any of the detunings for relevant lattice depths 



(Vl ~ 20E r ). Thereby, all couplings apart from the one 
being resonant is well in the dispersive regime and should 
only contribute with virtual processes. In other words, it 
is justified to introduce an effective 2x2 model for the 
<pi(t) and ip%(t) states alone; 
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L 12 u 





(11) 



showing perfect Rabi oscillations between the two states. 
The inverse Rabi frequency ? 12 is the effective coupling 
obtained when all the dispersive states have been adia- 
batically eliminated, i.e. in lowest order (in terms of cou- 
pling strengths divided by detunings) we have t s 12 = t^ 2 . 

Equation (ITTj) is the well-known Rabi model and is 
often considered for bosonic Josephson systems. It is 
known that interaction among the bosons (atoms) may 
greatly affect the Josephson dynamics in terms of self- 
trapping [l9| . Assuming lowest order in the coupling 
parameter and including onsite atom-atom interaction 
results in the effective model 
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where U s = U* 



cincl Up - Upppp 



(12) 

are the onsite effective 
scattering coefficients. Whenever the diagonal elements 
U 8 \(ff (t)\ 2 and U p \ip^(t)\ 2 are different, the dynamics is 
efficiently detuned and the population transfer is not per- 
fect. Thus, the interaction can been seen as inducing 
effective shifts of the onsite energies. This self-trapping 
effect is demonstrated in Fig. [5] displaying the popula- 
tions P a (t) = |</9 Q (£)| 2 ; in (a) the atom-atom interaction 
is zero while in (b) it is non-zero. We show the numerical 
solutions of the two-site six-state model and the solutions 
of the effective two-state model (fT2")) . Without interac- 
tion, the agreement between the two models is almost 
perfect and therefore the transfer is approximately 100 
%. The results disagree somewhat when the interaction 
is non-zero, even if the qualitative structure is the same. 
Within the six-state model we included all non-zero in- 
teraction terms Ua^s- The effective interaction strength 
U s Nq used in (b) is very large in order to demonstrate 
the self-trapping phenomena. In fact, we have in this 
example U s Nq ~ 4.8|tf|| implying the strongly interact- 
ing regime. For example, the critical coupling for en- 
tering the Mott insulator state with one atom per site 
U s Nq w 4|tf|| [20j - For a factor 10 smaller coupling, the 
six-state model predicts a transfer of 96.7 %, and for a 
factor 100 the transfer is 98.8 %. For typical experimen- 
tal parameters it is hence expected that the self-trapping 
effect will be vanishingly small, and in the following sub- 
section we will focus on a non-interacting gas. 



B. Full one dimensional dynamics 

We showed above that in the two-site system, an effec- 
tive two-level model qualitatively describes the dynam- 
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FIG. 3: Maximum population transfer P m ax(c) as a function 
of the dimensionless parameter c. The solid line is the result 
for one dimension and the dotted for two dimensions. The 
potential depth Vt = 20E r and U) = 40 x 2ty Hz. 



ics. In a TV-site system there are numerous of partic- 
ipating levels, even though the full Hamiltonian is six- 
diagonal. Adiabatically eliminating all far detuned levels 
might therefore induce non- vanishing shifts in the result- 
ing effective two-level model. Assuming a two-level model 
to be valid, we modify Eq. (TTTI) to allow for a detuning 
between the two levels 



E-2 — E\ 



ih 



(So - c) 



<p{(c,t) 
<f^(c,t) 



.Eg — Ei 
2 



(So 



(pl(c,t) 
f p 2 (c,t) 



(13) 

Thus, if the full TV-site Hamiltonian shows similar dy- 
namics as predicted by (|13p we conclude that the system 
Hamiltonian is effectively 2x2 block diagonal. More- 
over, for vanishing shifts g — and So — 1 so that 
c = 1 result in a vanishing effective detuning (remember 
that c was introduced in order to adjust for any virtual 
energy shifts). Denoting A(c) = (E 2 — Ei)(5 — c)/2, 
the maximum population transfer P m ax(c) to the p-band, 
provided ip s (c, t = 0) = 1 , is given by the Lorentzian 



Ac) 



A 2 (c)+ 5 2 



(14) 



having a width ~ g. In Fig. [3] (solid line), we test this hy- 
pothesis by calculating the maximum population transfer 
for the TV = 41 site model including all 67V states and a 
harmonic trap. For the numerics, we truncate the lattice 
such that the atomic density is approximately zero (due 
to the confining trap) at the edge of the lattice. The ini- 
tial state is taken as the ground state for the F = case. 
Pmax(c) m Fig- El shows a clear Lorentzian shape indi- 
cating that the dynamics could be effectively described 
by a two-level model. Corrections are manifested in, for 
example, the fact that P m ax(c) does not exactly reach 
1. The offset <5 (« 0.99998) from unity gives the shift 



in the onsite energies, while the width determines the 
effective coupling between the s- and the p-band. Fit- 
ting a Lorentzian to the data of Fig. [3J gives a coupling 
g « 0.99^2- The fact that g is only slightly shifted from 
t*2 and that approximates unity indicate minimal vir- 
tual shifts arising from off-resonant levels. 

This far, we have not included any coupling of the res- 
onances to the continuum, i.e. taking the imaginary part 
r„ of the Wannier-Stark energies into account. Utilizing 
a wave- packet propagation method [la ], we estimate T v 
by calculating the decay of population within the first two 
Brillouin zones as time progresses. This gives us a decay 
rate r = Y/% rs 0.1 s" 1 , giving a survival time (~ 10 s) 
which is considerably larger than the characteristic Rabi 
time (~ 0.1 s). Another aspect that we have not been 
investigated is the effects deriving from switching on/off 
of the lattice acceleration. Figure [5] shows that half a 
Rabi period is approximately given by T r w 27r7i/|i 1 2|, 
which we will use for the pulse duration in the expres- 
sion for F(t) in Eq. ([6|). For minimizing the loading time, 
T s should be chosen small. However, a too small switch 
on/off time T s may cause non-adiabatic excitations. We 
pick T s = 10 //s, giving F(t) a square pulse shape with 
rounded edges. The result of our numerical solution of 
the full one dimensional system including a trapping po- 
tential and a pulse shaped acceleration is presented in 
Fig-Hl The final population on the p-band is 95 %. The 
effect of the imaginary coupling to the continum is a drop 
of the population transfer by approximately 1 %. 
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FIG. 4: Populations Pi(t) (i = s, p, d) for the pulsed acceler- 
ation @ with T s = 10 /Lts and T r = 27rR/|t*f +1 |. The number 
of sites was taken as N = 41, well beyond the extent of the 
atomic density. The potential depth Vl = 20E r , c = 1, the 
loss rates T v /fi = 0.1 Hz (y = s, p, d), and cj = 40 x 2tt Hz, 
while the interactions are set to zero. The population transfer 
from the s- to the p-band is ~ 95.5 %. 



IV. LOADING IN A TWO DIMENSIONAL 
SYSTEM 

The method used for the one dimensional loading 
transforms directly to higher dimensions. In Fig. [31 the 
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dotted curve displays the maximum population trans- 
fer to the p-bands in two dimensions. The initial state 
is again the ground state for F — and the trap fre- 
quency lu — 40 x 2ir Hz. The shape is still Lorentzian, 
and the shift 5a f=a 0.99990. However, the effective cou- 
pling is roughly four times as large as in one dimension (a 
Lorentzian fit gives g 3.57ti2), which implies a speed- 
up of the p-band loading. 

One difference between the two cases is that in two 
dimensions there are two p- and three ci-bands. In addi- 
tion, atoms can hop in two directions, x and y. If an atom 
in a s-state hops to a p-state orbital in the ^-direction, 
due to symmetry it can only hop into a p-state orbital 
with a node in the z-direction. Similarly for hopping in 
the y-direction. Thus, neglecting all off-resonant transi- 
tions and atom-atom interaction one derives an effective 
three-state model 



lh di 



where we have labeled the three sites 1, 2, 3. By defining 
(p s (t) =(f1{t) and ip p ± (t) = (<p%(t)±<pl(t))/V2axLd using 
t sp we find an effective two-site model 
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Thus, comparable to the one dimensional loading, here 
the coupling has been increased by a factor However, 
Fig. [3] made clear that the effective coupling is in fact in- 
creased by a factor ~ 4 rather than y/2. This is assumed 
to derive from the fact that there are many more off- 
resonant levels and the corresponding virtual processes 
influence the effective coupling more drastically when all 
of them are eliminated. 

Since there are multiple p- and d-bands we define the 
total populations as 



P P (t) 
Pd(t) 
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P X (t)+Py(t), 
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band 



band 



(17) 



Po y (t), d — band. 



We will mainly be interested in the transfer efficiency and 
not in the actual quantum state loaded in the p-bands. 
This can be motivated from the fact that the prepared 
state is not a stationary one Q , and non-trivial dynamics 
continues also when the acceleration has been turned off. 
For example, in the experiment of Ref. Q it was shown 
how condensation appeared on the p-bands after loading 
into these bands. 

For loading of p-band atoms we again use the pulsed ac- 
celeration (|6|) but this time with c = 0.9999. We see from 
Fig-Elthat the corresponding two dimensional Lorentzian 
is shifted away from c = 1 more than the one dimensional 
Lorentzian. In the previous example we could keep c = 1 
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FIG. 5: Probabilities P a (t) and P p (t). The final transfer to the 
p-bands is 96.5 % in the non-interacting case (black curves) 
and 83.0 % for the interacting situation (pink curves). The 
parameters are the same as in Fig. [4j but the dimension- 
less parameter c = 0.9999, and the half Rabi oscillation time 
was picked as T r = 39 ms. The interaction strengths have 
been taken very high in order to achieve a self-trapping effect; 
(TooooiVo w 12\t av \ and U XXXX N = U ym N » 51|t sp |. Losses 
into the continum is the same as in Fig. 01 i.e. T v /h = 0.1 Hz 
(v = 0, x, y, xy, xO, Oy). 



and still attain a large population transfer. In two di- 
mensions, on the other hand, utilizing c = 0.9999 in- 
stead of c = 1 increases the efficiency by several percent. 
Similarly, since the virtual processes are more impor- 
tant in this two dimensional situation we do not chose 
T r = \/2t sp , but find T r numerically. The results for 
the full system evolution including the trap are shown in 
Fig. [5] as black curves. Losses r„ are taken the same as 
for the one dimensional case. During the pulse switch 
off, the p-band population decreases with half a percent 
indicating non-adiabatic effects. 

Our numerical simulations allow us to extract other 
quantities, like the densities n x (i x ,i y ) and n y (i x ,i y ) 
which we display in Fig. [6j Interestingly, the two den- 
sities are not identical and the y-density is slightly (5 
%) less populated. The origin of this asymmetry is not 
clear. Nevertheless, we note that since the two p-bands 
are degenerate and an imbalance of population between 
the two flavors does not change the energy of the system, 
i.e. the results are very sensitive to any small fluctua- 
tions. Atom-atom interaction might stabilize this, but as 
we are not focusing on the actual quantum state on the 
p-bands in this work, we do not examine this further. 

In the interacting case, we do not solve the full system 
which would contain more than 1000 different interac- 
tion terms. Instead we use the fact that mainly the s- 
and p-bands are populated and therefore only include the 
interaction terms 



U, 



0000, 



Ur 



yyyyi 



(18) 



corresponding to scattering processes that preserve the 
atomic orbital structure. Interactions that do not sus- 
tain orbital order, like Uqq xx , Uoo yy , and U xxyy , are all 
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smaller than the above. We find that moderate interac- 
tion can actually increase the population transfer. How- 
ever, for extremely strong interaction, like in the one 
dimensional situation, a self-trapping effect sets in. In 
two dimensions, the system is less sensitive to interac- 
tion and therefore self-trapping becomes evident only for 
extremely strong atom-atom interactions. This is dis- 
played in Fig. [5] as the pink curves. In this example, 
the strength of the interactions, U a pj5, are one order 
of magnitude larger than the tunneling coefficient t sp . 
The transition where interaction induced self-trapping 
becomes important is found to be relatively sharp. For 
example in Fig. for 20 % smaller interaction strengths 
the population transfer is increased to around 96 %. 




lated but their mutual dynamics is complex. In this ir- 
regular regime, at some occasions the p-band population 
can become very high while in the next moment decrease 
rapidly. At Vl ~ 18E r , the dynamics become more Rabi- 
like and from there on we encounter an increased loading 
efficiency with growing Vj,'s. Experimentally one should 
pick a large Vl in order to have a stable loading, but too 
large Vl will imply long loading times causing heating 
and other incoherent processes to become important. 



0.99 



Oh 




14 16 18 20 22 24 
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r 



FIG. 7: The maximum population transfer P max from the s- 
to the p-bands as a function of the potential depth Vl ■ The 
interaction between the atoms is assumed zero, and the decay 
rate is taken the same for all E r and its value is as in Fig. [4] 




FIG. 6: Normalized atomic densities n x (i x ,i y )/No (a) and 
n y (i x ,i y )/No (b) on the p-band. As expected, the densities 
are mainly localized around the trap center. 



We finish this section by investigating the transfer effi- 
ciency in terms of the potential depth. We consider non- 
interacting atoms and calculate the maximum population 
Pmax on the p-bands for various Vl- The results are pre- 
sented in Fig. [7] At smaller depths, we have noticed 
that the population transfer between the bands are no 
longer displaying a Rabi-like oscillation structure. Even 
at Vl = 10E r , mainly the s- and the p-bands are popu- 



V. CONCLUSIONS 

We have shown that by properly accelerating optical 
lattices, bosonic atoms can be transfered from the low- 
est s-band to the first excited p-bands. The idea re- 
lies on creating a resonant coupling between the s- and 
the p-bands in the lattice. Since other transitions are 
off-resonant, the population mainly resides within these 
bands. With an efficiency of more than 95 % population 
transfer, the scheme seems as a good alternative com- 
pared to earlier methods like Raman pulses [13] , ramped 
super-lattices [1, H|, or condensation into moving lat- 
tices [7j. 

The influence of atom-atom interactions was studied, 
and it was found that for strongly interacting gases self- 
trapping may deteriorate the p-band loading. This effect 
is considerably reduced in higher dimensions. Even in one 
dimension, the interaction must be very large in order to 
become influential. Experimentally, the interaction can 
be made weaker by the use of Feshbach resonances [2l[ 
or opening up the trap, e.g. we consider a trapping fre- 
quency ui ~ 250 Hz and it could in principle be decreased 
by one order of magnitude [7]. Tunneling into the con- 
tinuum was as well estimated to have a minor influence 
on the loading. In particular, for such deep lattices as 
20E r which we consider, the decay rate was found as 
0.1 s _1 which is an order of magnitude smaller than the 
characteristic loading rates. 
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Throughout, we have tried to employ experimentally 
relevant parameters. For example, we considered 87 Rb 
atoms and a lattice depth Vl — 20E r . This resulted in a 
loading-time of roughly 40 ms in two dimensions, which is 
within the coherence and heating times of a BEC [8j . For 
the given parameters, the acceleration a = Fn/m « 500 
m/s 2 which is easily achieved experimentally [7j. The 
ability of the method for other lattice depths was as well 
discussed. Despite the fact that we consider experimen- 
tally relevant parameters and well developed techniques, 
there are possible sources that may lower the efficiency. 
For example, the method has a non-adiabatic character 
and it therefore relies on correctly adjusting the acceler- 
ation in terms of strength and duration. It follows from 
Fig. [3] that this can become a serious issue in one dimen- 
sion where the population transfer is extra sensitive to 
the pulse strength. In two dimensions, this is less cru- 
cial, but one would probably need a control of the pulse 
strength within less than one percent. Another impor- 
tant aspect is the specific strengths of the tunneling co- 
efficients obtained numerically from Wannier functions. 
The actual form of the Wannier functions, and especially 
their tails, greatly affects the values of the overlap inte- 
grals. Experimentally, the system might not occupy a 
true stationary state at t = and such corrections might 
affect the various tfJ 3 . It is also known that atom-atom 
interaction modifies the Wannier functions [22j, which 



will propagate into the tunneling coefficients. 

We considered a tight-binding approach and expanded 
the atomic field in terms of the three lowest band's Wan- 
nier functions. For the s-band, the tight-binding approx- 
imation is typically justified for Vl > 3E r [23J. For the 
higher excited bands, the validity of the tight-binding 
approximation is more stringent. This, however, is not 
assumed to cause any problems in our analysis since the 
higher bands are off-resonant beyond nearest neighbor- 
ing sites, and even if they couple to these sites the effect 
is greatly suppressed by the large detuning. Therefore, 
by carefully controlling the experimental parameters, the 
present loading scheme seems to be experimentally feasi- 
ble with current technology. Furthermore, we note espe- 
cially two most recent experimental works [24| in which 
resonant tunneling between neighboring sites has been 
studied. Contrary to our work, Ref. [24| considers a 
strongly interacting system in the insulating regime, but 
it nevertheless demonstrate coherent population transfer 
between neighboring sites in tilted lattices. 
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